Development of new materials for electrothermal metals using data driven and machine learning

After adopting a combined approach of data-driven methods and machine learning, the prediction of material performance and the optimization of composition design can significantly reduce the development time of materials at a lower cost. In this research, we employed four machine learning algorithms, including linear regression, ridge regression, support vector regression, and backpropagation neural networks, to develop predictive models for the electrical performance data of titanium alloys. Our focus was on two key objectives: resistivity and the temperature coefficient of resistance (TCR). Subsequently, leveraging the results of feature selection, we conducted an analysis to discern the impact of alloying elements on these two electrical properties.The prediction results indicate that for the resistivity data prediction task, the radial basis function kernel-based support vector machine model performs the best, with a correlation coefficient above 0.995 and a percentage error within 2%, demonstrating high predictive capability. For the TCR data prediction task, the best-performing model is a backpropagation neural network with two hidden layers, also with a correlation coefficient above 0.995 and a percentage error within 3%, demonstrating good generalization ability. The feature selection results using random forest and Xgboost indicate that Al and Zr have a significant positive effect on resistivity, while Al, Zr, and V have a significant negative effect on TCR. The conclusion of the composition optimization design suggests that to achieve both high resistivity and TCR, it is recommended to set the Al content in the range of 1.5% to 2% and the Zr content in the range of 2.5% to 3%.


Introduction
In the realm of new material development, the challenge of achieving specific target performance has spurred a revolutionary transformation in material research, rendering traditional trial-and-error methods insufficient.These conventional approaches rely on limited experimental data to discretely adjust material compositions and process parameters in the quest for optimal material performance.Nevertheless, it is evident that this research and development methodology is time-consuming and costly.In an effort to expedite material development and introduce more efficient research approaches, the United States initiated the Materials Genome Initiative (MGI) in 2011 [1], with the aim of significantly reducing the material development cycle.
In recent years, data-driven and machine learning methods have gradually emerged as prominent players in addressing real-world engineering challenges [2][3][4][5][6][7][8][9][10].The application of artificial intelligence, machine learning, and deep learning technologies in the field of materials science has garnered significant attention, as these methods are widely adopted for material discovery and design [11].Logan et al. [12] emphasized that the process of using machine learning methods in materials informatics comprises three key elements: material data, material descriptors, and machine learning algorithms suitable for prediction.Davoodi et al. [13] employed multiple machine learning models to accurately predict the hydrogen (H2) absorption percentage in porous carbon media (PCM), offering critical insights for efficient H2 storage.Bruno et al. [14] utilized a machine learning kernel regression model to predict the electronic properties and elastic performance of materials, along with proposing an Exhaustive Enumeration algorithm for material reverse design.The majority of the methods mentioned above have effectively utilized suitable machine learning algorithms, coupled with selective feature engineering, to achieve satisfactory prediction and design outcomes.
The objective of this study is to develop electrically conductive metallic materials with high resistivity and temperature coefficient of resistance (TCR).In the domain of alloy electrical performance, while research on resistivity is extensive, investigations into the temperature coefficient of resistance (TCR) of alloys remain relatively limited.However, both resistivity and TCR are critical indicators of alloy electrical performance.Faced with the vast parameter space of multicomponent alloys, traditional experimental methods are practically incapable of encompassing all conceivable alloy combinations.Consequently, this study introduces datadriven and machine learning approaches to address the challenges of compositional optimization design, thereby accelerating the development of new materials.
The remaining sections of the introduction are as follows: Section 2 will introduce various machine learning algorithm models, encompassing predictive modeling and feature selection.Section 3 will compare a series of experimental results to determine the optimal predictive model and discuss and analyze compositional optimization.Lastly, Section 4 will summarize the work and propose future research directions.

Method
This paper presents a three-step data-driven and machine learning-based workflow for predicting the electrical performance of electrically conductive metal materials.As Fig 1 illustrated, in the first step, the electrical performance data was obtained and subjected to data cleaning, followed by the partitioning of the original dataset into a training set and a test set.Subsequently, the raw data was normalized and standardized to mitigate the impact of varying data scales on the prediction results.In this study, we employed K-fold cross-validation as the method for data partitioning and model training, which has been proved as a well-established strategy [15][16][17].K-fold cross-validation was utilized to reduce the influence of different data partitions on the training process.In this study, we set K to 5 for dataset splitting.The division ratio of 4:1 provides us with approximately 80% of the data for training, leaving 20% of the data for testing.Within the training data, a further 7:3 split is used to create a validation set, which is employed for assessing the model's performance and evaluating its ability to generalize to unseen data.Separating the validation set from the training set allows for a more robust estimation of the model's performance in practical applications.This balanced approach ensures that while a sufficient amount of data is retained for model performance evaluation, there is also ample training data to support the model's learning, thereby safeguarding its generalization capabilities.Machine learning methods exhibited varying sensitivities to material data within different ranges [18][19][20][21].Hence, it was crucial to select proper algorithms based on specific material data samples and evaluate their performance using appropriate performance evaluation metrics.In this study, we evaluated the models using four metrics: root mean square error (RMSE), mean absolute error (MAE), mean absolute percentage error (MAPE), and coefficient of determination (R 2 ).Among these, R 2 serves as the principal performance metric to identify the most effective predictive model for subsequent analysis.
The second step involved feature importance selection.In this study, we utilized two different feature selection models and two evaluation methods to rank the importance of features.The intersection of the feature subsets obtained from both methods is selected as the final feature subset [22].Since the feature data in this study comprised elemental composition data, the results of the feature importance selection represented the degree of influence that different elements have on the electrical performance of alloys.As our target performance included resistivity and TCR, separate feature selection was performed for each of these electrical performance indicators to identify the element subsets that have the greatest impacts on each of them.
The third step utilized the outcomes of the previous two steps for compositional optimization analysis.Due to the "trade-off" relationship between resistivity and TCR, achieving simultaneous positive gains in both electrical performance indicators cannot rely solely on blindly increasing or decreasing the content of certain elements based on feature selection results.In this study, by adjusting the element compositions obtained through feature selection, novel electrical performance data was calculated using the predictive model.The trends in resistivity and TCR changes were then comprehensively analyzed within the same range to explore the design principles for optimizing these two electrical performance indicators simultaneously.The details of each step will be discussed in the following subsections.involved four main processes.Firstly, various types of alloy electrical performance data were collected through literature review, wherein the alloys were initially selected based on the values of resistivity and TCR.In this study, titanium alloys were chosen as the initial research objects.Next, programming was employed to generate 729 virtual sample points with titanium as the primary constituent element.These 729 titanium-based compositions included specific elements such as Al, Si, Zr, Mo, V, Sn, Nb, Mn, and Ti, totaling nine elements.Subsequently, the JMatPro high throughput modules (named API), which allows batch calculation, was utilized to obtain the resistivity data at 29 temperature points ranging from 25˚C to 300˚C for each virtual composition sample.A total of 21,141 resistivity data points were obtained for titanium alloy, constituting the resistivity dataset used in the following experiments.Lastly, by utilizing the calculation relationship between resistivity and TCR, a programming approach was employed to calculate the TCR data for the 729 compositions, forming the TCR dataset used in the following experiments.

Calculation of TCR
With 25˚C as the reference temperature, the relation between the resistance R t of the alloy at t ˚C and the temperature can be expressed by the following formula: In Formula (1), R 25 is the resistance value at 25˚C, and the unit is O; α, β and γ are the resistance temperatures of 1st, 2nd and 3rd power respectively.Since the resistance value and temperature of titanium alloy in this paper are approximately linear in the range of 25˚C ~300˚C, the relationship between resistance and temperature can be simplified by Eq (1) as follows: From the relationship between resistance and resistivity, it can be seen that the change of resistance value of the same material at different temperatures essentially comes from the change of resistivity, so Formula (2) can be further simplified as: In Formula (3), ρ t and ρ 25 are the resistivity of titanium alloy at t˚Cand 25˚C respectively, and the unit is O�m.As shown in Fig 2, t-25 is the x-axis and ρ t is the y-axis for linear fitting.The fitted intercept is ρ 25 and the slope is ρ 25 α, so the resistance temperature coefficient α can be calculated by dividing the slope by the intercept.

Data normalization and standardization
Data normalization and standardization can also keep the original distribution of the data while eliminating the influence of the large difference in the order of magnitude of different features on the prediction model.Therefore, this study will normalize or standardize the original data before regression prediction and feature screening.
In normalization processing, each feature is transformed by the following formula: Where: x ¼ ðx 1 ; x 2 ; � � � ; x m Þ represents the set consisting of the values of a feature of m samples; x i 2x; max(x) and min(x) are the maximum and minimum values in the feature set respectively.x * i indicates that the eigenvalue of x i in any value range is transformed into the value in the interval [0,1] through calculation, that is, the eigenvalue after normalization.
In normalization, each feature is transformed by the following formula: Where, mean(x) and std(x) are the characteristic mean and standard deviation respectively.The mean value of the data is 0 and the standard deviation is 1.The reason for the above data processing is that after comparing the two methods, it is found that the prediction effect after standardized processing is better.

Machine learning regression model
In predictive tasks, machine learning can be divided into classification and regression [23][24][25].The prediction of electrical properties in this paper is a regression task, wherein, four commonly-used machine learning models were employed: linear regression, ridge regression, support vector regression and BP neural network.Among them, BP neural network showed strong nonlinear mapping ability and was widely used in classification, fitting, diagnosis, prediction and other fields [26].
Linear regression.Linear regression model plays an important role in the field of machine learning, and its advantages are simple and easy to model.Its basic form is: Let the given data set number of samples, d is the feature dimension, y i 2R.Linear regression tries to learn: The linear model has a certain interpretability.It can be seen from Eq (6) that ω i represents the weight of each feature in the prediction, and b is the target value when all features are zero.When using the least square method to solve ω and b, that is: The optimal solution of ω and b can be obtained from Eq (7).Ridge regression.Ridge regression is generally used in cases where the ratio of sample characteristic dimension to sample number is large, and it is a supplementary method to linear regression.The loss function of ridge regression model is: L2 norm regularization is used in the expression, and the regularization parameter λ>0.
Ridge regression can alleviate the overfitting problem by introducing regularization terms.In addition, from A mathematical point of view, Ridge regression can better solve the problem that linear regression can not get more stable ω or can not be solved.
Support vector regression.The principle of support vector machine [27,28] is to find a hyperplane to divide different classes of samples.The objective function of the SVR problem is as follows: Where: ω = (ω 1 ,ω 2 ,� � �,ω d ) is the normal vector dividing the hyperplane; C is the regularization constant; l ε is the loss function.
SVR can be expressed by the formula with kernel function as: Where, κ(x,x i ) is the kernel function, ^ai and a i are Lagrangian multipliers of the ith sample, ^ai �0, a i � 0. The SVR cores used in this study are Gaussian radial basis (RBF kernel) and S-type kernel (Sigmod).
BP neural network.BP neural network is the extension of perceptron, the main characteristics of this network are the forward transmission of working signal and the reverse propagation of error signal.The internal network layer structure can be divided into three parts, namely, the input layer, the hidden layer and the output layer, as shown in Fig 3 .Each layer of the network contains a number of neurons, and the layers are connected in a fully connected way, and the output of the upper layer is the input of the next layer, so as to realize the feature mapping and mining the inherent law of the data.
Suppose that the number of neurons in layer l is n, then the output vector of layer l is α l 2R n×l .The output of layer l can be expressed by the formula: Where: z l 2R n×l is the input vector of layer l; W l 2 R n lÀ 1 �n l is the weight matrix from layer l-1 to layer l; b l 2R n×1 is the bias vector of the l layer; σ(�) is the activation function.Common activation functions include Relu, Tanh, Sig-moid and so on.
Evaluation index.The root mean square error (RMSE), mean absolute error (MAE), mean absolute percentage error (MAPE) and goodness of Fit (R 2 ) were used to evaluate the performance of the prediction model.
RMSE is calculated as follows: RMSE ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi 1 m MAE is calculated as follows: MAPE is calculated as follows: The maximum value of R 2 is 1.The closer the value is to 1, the better the fitting degree is.The calculation method is as follows: Among the above types, y i and ŷi are the real and predicted values of the i sample respectively, with the mean value ŷ ¼ 1.
Grid parameter optimization.In the realm of machine learning, models typically comprise numerous hyperparameters that require careful tuning to achieve optimal performance.These hyperparameters encompass factors such as learning rate, the number of trees, depth, and regularization strength.The selection of the right combination of hyperparameters is pivotal for attaining high-performance models.
The fundamental principle of grid search involves exploring a predefined grid of hyperparameter value combinations.Initially, a grid of hyperparameter combinations is constructed based on the hyperparameters and their candidate values.It is essential to exercise caution when choosing hyperparameters and candidate values to prevent the search space from becoming excessively large, consuming excessive computational resources.Subsequently, for each hyperparameter combination, models are trained on the training data, and their performance is evaluated using techniques such as cross-validation.Finally, based on the results of performance metrics, the hyperparameter combination exhibiting the best performance is selected.This combination typically exhibits the lowest error or the highest accuracy.
This study necessitates the use of three models for grid search hyperparameter optimization: BPNN, Support Vector Machine, and Ridge Regression.The grid parameter settings for these three models are as follows:

Feature selection
In traditional regression analysis models, the more features means the less accuracy the predictions become.Therefore, feature selection is necessary to improve results and optimization.
Feature engineering is a crucial step in machine learning [29][30][31], and different feature choices can significantly impact the prediction outcomes.In the context of titanium alloy electrical performance prediction, when the prediction model utilized elemental features to forecast target electrical performance, each elemental feature became a factor influencing the prediction results.The task of feature selection is to identify highly correlated features with the target performance while filtering out irrelevant features, thereby improving prediction effectiveness and reducing computational complexity.The commonly used feature selection methods mainly include three types: filter-based, embedded, and wrapper-based approaches [32].In this study, a wrapper-based feature selection method was employed, which combined feature subsets with machine learning models and utilized the predictive performance of the models as the criterion for selecting the feature subsets.Since the feature selection process involves multiple iterations of training the learner, wrapper-based feature selection entails significant computational complexity and time consumption.This paper employed two wrapper-based feature selection algorithms, namely random forest and Xgboost, with feature importance and SHAP (SHapley Additive exPlanations) values used to illustrate the feature selection results.
Random forest feature selection.Random forest adopts random feature selection method to reduce the overfitting risk and increase the generalization ability of the model by randomly selecting different feature subsets during the construction of each decision tree.
The random forest defines the importance measure of X feature.The steps to calculate the importance of a feature are as follows: 1.The prediction Error rate of random forest for samples outside the bag is called Out-Of-Bag error (OOB).For decision tree T i in random forest, the number of classification errors E i on the data outside the bag is calculated.
2. The value of X is randomly disturbed in the OOB data of the decision tree, and the number of classification errors E X i is recalculated.3. Let i = 1,2,. ..,n, repeat the above two steps, where n is the number of decision trees contained in the random forest.
4. The importance of feature X is defined as: The basis of this definition is that if the OOB error of the model is significantly increased after the addition of noise to a feature, it indicates that the feature has a greater impact on the prediction result, and thus has a higher importance.
Xgboost feature selection.Xgboost uses incremental training and fine-tuning of feature split points, so that variable selection and weight parameter adjustment can be optimized in continuous iteration to improve model performance.The core idea of the Xgboost algorithm is to generate trees by constantly splitting features, and for each tree generated, a new function is generated to fit the residual of the last prediction.When we need to predict a sample, each sample will fall on the corresponding leaf node in each tree, and will correspond to a score, and finally add the scores corresponding to each tree to get the predicted value of the sample.
The objective function of Xgboost algorithm is: In Eq (18), l is the loss function, O(f t ) is the regular term, and const is the constant term.SHAP.SHAP is a widely recognized attribution method in the field of explainable machine learning, derived from the Shapley value in cooperative game theory [33].The key to interpreting machine learning models using SHAP analysis is to calculate the Shapley value for each feature of every instance.In practical applications, the Shapley value is approximated, and SHAP is an optimization algorithm for estimating Shapley values.
This method possesses three desirable properties of Shapley value theory: local accuracy, missingness, and consistency.Local accuracy refers to the sum of feature attributions being equal to the difference between the expected output and the actual output of the explanatory model, enabling the explanatory model to capture the discrepancy between the desired output and the observed output of a given instance [34].Missingness implies that the Shapley value for a feature that is absent is considered as 0, indicating no contribution to the prediction.
Consistency means that if a change in the model leads to an increase or unchanged marginal contribution of a particular feature, the Shapley value for that feature will not decrease.In addition to these features, the greatest advantage of SHAP lies in its independence from the predictive model [35], allowing for the theoretical interpretation of any machine learning model.

Results and analysis
The target performance data and feature data for the experiment were obtained by fitting and calculating 21141 titanium alloy resistivity samples using Eq (3), resulting in 729 sets of TCR (Temperature Coefficient of Resistance) data.The resistivity data at 25˚C for the 729 sets of titanium alloys, as well as the composition proportion data containing nine elements, constituted the experimental feature data.The target performance data and feature data together formed the experimental dataset, which is summarized in Table 1.And the statistical analysis results for Resistivity and TCR are individually presented in (a) and (b) of Fig 4.
In Table 1, "Count" represents the total number of data points for each data type, which is 729 sets in this case."Mean" denotes the average value of each data type, while "Std" represents the standard deviation."Min" and "Max" correspond to the minimum and maximum values, respectively.For elemental compositions, they also indicate the compositional range of each element.It is evident that the titanium (Ti) content, as the base element for titanium alloys, ranges from 66.28% to 98.53%, making it the most abundant element among all.In Fig 4(A), the histogram represents the distribution of resistivity values, with the vertical axis indicating the frequency of resistivity values.The purple curve represents the kernel density estimation curve.From the frequency histogram, it is evident that the primary range of resistivity values falls between 0.7 and 1.6.The kernel density estimation curve exhibits three peaks, corresponding to intervals with a higher concentration of data points, namely around 0.9, 1.After normalizing or standardizing the data in Table 1 through preprocessing operations, subsequent prediction modeling, feature selection, and optimization range analysis were performed using the resistivity and TCR as target performance data.

Predicted results of resistivity model
Table 2 presents the predicted results of resistivity data based on five machine learning models, and the values of four evaluation metrics are based on the average level of K-fold cross- validation.The BPNN (Backpropagation Neural Network) model and three SVR (Support Vector Regression) models with different kernel functions are displayed with their best results after parameters tuning.
From the results shown in Table 2, it can be observed that the Rbf (Radial Basis Function) kernel SVR model performs the best among all models, with superior results in all four evaluation metrics, and the coefficient of determination R2 exceeding 0.99.
The parameters of the Rbf-SVR model for resistivity prediction are as follows: The Rbf kernel is chosen for SVR, with a regularization coefficient C of 150, a kernel coefficient gamma of 0.1, and ε set to 0.01.In Fig 5(A) to 5(E), the regression performance of five machine learning models on the validation and test datasets is presented.It is evident that the SVR model based on the Rbf kernel demonstrates the best regression performance, with minimal discrepancies between the predicted and actual resistivity values, resulting in an overall excellent predictive performance.In conclusion, the Rbf-SVR model is selected as the predictive model for resistivity, the target property in this study.

Prediction results of TCR model
Table 3 presents the prediction results of TCR data based on five machine learning models, and the values of the four evaluation metrics are still based on the average level of K-fold crossvalidation.Similarly, the BPNN model and three SVR models with different kernel functions are displayed with their best results obtained during parameter tuning.
From the results shown in Table 3, it can be observed that the BPNN model achieves the best prediction performance among all models, with superior results in all four evaluation metrics, especially with a coefficient of determination R2 exceeding 0.99.The Rbf kernel SVR model follows closely in terms of prediction performance, but the errors, as indicated by MAE and MAPE, are approximately twice that of the BPNN model.The linear regression, ridge regression, and Sigmod kernel SVR models exhibit similar prediction performance, however, their performance is significantly inferior compared to the BPNN model.
The parameters of the BPNN model for TCR prediction are as follows: The neural network has two hidden layers, each with 36 nodes, and the activation function is set as Relu.The weight optimizer is set to stochastic gradient descent, with a regularization parameter of 0.0001, and the maximum number of iterations is set to 10,000.
In  regression performance, with minimal discrepancies between the predicted and actual TCR values, resulting in an overall excellent predictive performance.In conclusion, the BPNN model is selected as the predictive model for TCR, the target property in this study.

Resistivity feature selection
Since the Ti element is the base element and its content often varies with the content of other trace elements.As a result, Ti is not involved in the feature selection process.The Si and Mn elements are trace elements with fixed composition, so these two elements are also excluded from the feature selection process.The impact of different elements on the resistivity target performance was investigated using two different machine learning models, Random Forest and Xgboost, and the results are shown in Fig 7.  It is evident that the content of the Al element is highly significant for the model.Samples with higher Al element content, represented by red points with high SHAP values (greater than 0), exert a positive influence.Conversely, for the blue-shaded portion where feature values are smaller, SHAP values are negative (less than 0), indicating a negative influence.From a horizontal perspective, the distribution of samples for the Al element feature is more scattered, signifying a greater impact of this feature.Additionally, for features like Mo element, most data points are dispersed around SHAP = 0, suggesting that the Mo element has a minimal impact within the normal range.
Through the SHAP plots, it becomes apparent that both Al and Zr have a similar impact on resistivity.Specifically, within a certain range, higher concentrations of Al and Zr lead to increased resistivity.This insight provides valuable guidance for subsequent compositional optimization analysis.
TCR feature selection.Similarly, Ti, Si, and Mn are not involved in the feature selection process for TCR.The impact of different elements on the TCR target performance was investigated using the Random Forest and Xgboost models, and the results are shown in Fig 8.It can be observed that the impacts of Al, Zr, and V on TCR are opposite to that on resistivity.Within a certain range, increasing the content of Al, Zr, and V leads to decrease in TCR, which is consistent with the reciprocal relationship between resistivity and TCR during compositional design.

Composition optimization analysis
Based on the feature selection results for resistivity and TCR, the variable elements chosen for compositional optimization are Al and Zr.When adjusting the content of Al and Zr, the content of other elements is set to the average values, as listed in Table 1.The newly generated virtual samples are used as the test set, and the resistivity is predicted using the SVR model with the Rbf kernel and the TCR is predicted using the BPNN model with the optimal parameters.The trends of resistivity and TCR variations are then analyzed in the same coordinate system.
Effect of Al on electrical properties.The average values for Si, Zr, Mo, V, Sn, Nb, and Mn are set to 7.50%, 3.31%, 0.86%, 1.50%, 1.25%, 6.63%, and 1.40% respectively.According to Table 1, the composition range of Al is set to 0-5% with a step size of 0.1.The content of Ti changes with the content of Al.In total, 51 virtual samples are generated.
These 51 virtual samples are the test set for prediction using the Rbf-SVR and BPNN models.The resulting curves of resistivity and TCR prediction data are shown in Fig 9 .The red curve represents the resistivity variation caused by the change in Al content from 0% to 5%, while the green curve represents the TCR variation.It can be observed that as resistivity increases, TCR decreases significantly.Therefore, when designing the target electrical performance of titanium alloys, it is important to consider both resistivity and TCR variations.The yellow dashed line in Fig 9 .corresponds to an Al content of 3.8%.The intersection of the dashed line with the two curves represents the values of both electrical properties.It can be seen that at this point, the resistivity performance is good, but the TCR is low.To improve the TCR performance of the titanium alloy, the dashed line needs to shift to the left, indicating a decrease in Al content.Based on this optimization approach, it is evident that the intersection points of the two curves correspond to satisfactory levels of both resistivity and TCR.In this study, both resistivity and TCR are targeted for high performance.Therefore, the range of Al content can be considered within the vicinity of the intersection points shown in Fig 9, approximately around 1.5% to 2%.
Effect of Zr on electrical properties.The average values for Al, Si, Mo, V, Sn, Nb, and Mn are set to 2.56%, 7.50%, 0.86%, 1.50%, 1.25%, 6.63%, and 1.40% respectively.According to Table 1, the composition range of Zr is set to 0-6.6% with a step size of 0.1.The content of Ti changes with the content of Zr.In total, 67 virtual samples are generated.
These 67 virtual samples are the test set for prediction using the Rbf-SVR and BPNN models.The resulting curves of resistivity and TCR prediction data are shown in If the target performance leans towards either resistivity or TCR, this optimization approach can still be used for design purposes.

Fig 1 .
Fig 1. Overview of the modeling process in simplified form.After the optimal model is obtained, the variation trend of resistivity and TCR can be analyzed by adjusting the content of important elements based on the feature screening results.https://doi.org/10.1371/journal.pone.0297943.g001

2 , and 1 . 4 .
In Fig 4(B), the histogram depicts the distribution of TCR (Temperature Coefficient of Resistance) values, with the vertical axis showing the frequency of TCR values.The purple curve represents the kernel density estimation curve.The frequency histogram reveals that the main range of TCR values lies between 500 and 3000.The kernel density estimation curve indicates a gradual decrease in the number of samples with high TCR values.This reduction can be attributed to the fact that, in existing alloy materials, obtaining high TCR values is relatively more challenging compared to high resistivity values.

Fig 5 .
Fig 5. Cross-plots of simulated and predicted resistivity values for five machine learning models evaluated on training and testing subsets.Panels (a) through (e) correspond to the BPNN model, Rbf-SVR model, Sigmod-SVR model, linear regression model, and ridge regression model, respectively.https://doi.org/10.1371/journal.pone.0297943.g005 Fig 6.Cross-plots of simulated and predicted TCR values for five machine learning models evaluated on training and testing subsets.Panels (a) through (e) correspond to the BPNN model, Rbf-SVR model, Sigmod-SVR model, linear regression model, and ridge regression model, respectively.https://doi.org/10.1371/journal.pone.0297943.g006

Fig 7 (
A) and 7(B) present the feature selection results based on feature importance.From Fig 6(A), it can be observed that the top three elements ranked by feature importance using the Random Forest model are Al, Zr, and Nb.On the other hand, Fig 7(B) shows that the top three elements based on feature importance using the Xgboost model are Al, Zr, and Sn.Therefore, the feature selection for resistivity includes the common elements of feature importance from both machine learning methods: Al and Zr.In Fig 7(C) and 7(D) depict feature contribution representations based on SHAP values.In comparison to the feature importance method, which directly provides the importance ratio, the SHAP approach also reveals the polarity of the impact.The SHAP visualization results in (c) and (d) show close similarities.In the SHAP plots, each row represents a feature, with the x-axis indicating the SHAP values.The features are sorted based on the mean absolute value of

Fig 7 .
Fig 7. Characteristic screening results of resistivity.(a) and (b) represent the resistivity feature importance results based on random forest and Xgboost machine learning methods, respectively; (c) and (d) respectively represent the results of the two methods on SHAP.https://doi.org/10.1371/journal.pone.0297943.g007

Fig 8 (Fig 8 .
Fig 8. Characteristic screening results of CTR.(a) and (b) represent the TCR feature importance results based on random forest and Xgboost machine learning methods; (c) and (d) represent the results of the two methods on SHAP.https://doi.org/10.1371/journal.pone.0297943.g008

Fig 9 .
Fig 9. Effect of Al content change on electrical properties of titanium alloy.https://doi.org/10.1371/journal.pone.0297943.g009 Fig 10.The red curve represents the resistivity variation caused by the change in Zr content from 0% to 6.6%, while the green curve represents the TCR variation.The yellow dashed line in Fig 10 corresponds to a Zr content of 4.8%, where the resistivity performance is good and the TCR is low.Therefore, when designing the Zr content for the titanium alloy, the range can be considered around the intersection points shown in Fig 10, approximately around 2.5% to 3%.

Conclusions 1 .
The support vector machine model with the radial basis function kernel, Rbf-SVR, achieved the best performance in predicting the resistivity of titanium alloys, with a correlation coefficient of 0.995.The absolute percentage error between the true values and predicted values of the test samples is within 2%.The backpropagation neural network model with two hidden layers, BPNN, performed the best in predicting the TCR of titanium alloys, with a correlation coefficient of 0.996.The absolute percentage error between the true values and